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Abstract 

Recent experimental results for the ratio of the branching fractions of B — > D^tD t and 
r-j B —> D^'fiPp decays came as a surprise and lead to a discussion of possibility to constraining 

Q_i New Physics through these modes. Here we focus on B(B — > DtD t )/B(B — > Duv^) and 

argue that the result is consistent with the Standard Model within 2er, and that the test 
of compatibility of this ratio with the Standard Model can be done experimentally with a 
minimal theory input. We also show that these two decay channels can provide us with quite 
good constraints of the New Physics couplings. 

> 
o 

1. Introduction: Recent experimental result by BaBar pQ 

i fl <°>=f§5lSr a440±a058±0 - 042 - « 

seems to indicate a disagreement with the Standard Model (SM) prediction. In addition to the 
above ratio, the experimenters also measured the corresponding decays to the final vector meson. 
Since the latter involve more form factors and can also be experimentally more challenging due to 
the need of discerning the soft pion events in B — > D*(— > Dir)lv from those in B — > D**(— > Dir)lv, 
we prefer to focus on B — > DID decays, even if the issue of properly handling the soft photon 
emission in this decay still remains to be solved [2 J . [^] 

To begin with, let us write the relevant hadronic matrix element in the form in which it is 
usually done in QCD, namely, 

(D(p')\c^b\B(p)) = [p, +p',- ml ~ ml fr) F + (g 2 ) + m| " m * g^ ( g 2 ) , (2) 



1 Laboratoire de Physique Theorique est une unite mixte de recherche du CNRS, UMR 8627. 
2 £ labels the lepton flavor. In practice, I = fi refers to the combined B — > Dev e and B — > D^iv^, as both leptons 
can be considered as massless with respect to the heavy mesons involved in the process. 
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that leads to the differential decay rate, 
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From the above expressions it is obvious that for the massless lepton in the final state the scalar 
form factor F^q 2 ) does not contribute to the differential branching fraction. In the case of r- 
lepton, instead, the last term becomes more important and one can question whether or not the 
coupling to a scalar non-SM particle can be probed through this decay. Usual assumption is that 
a charged Higgs boson might give a non-zero contribution in the b — > cH~ — > cJLvi transition 0] , 
a scenario that was recently challenged by the experimental results . A couple of alternatives 
have already been proposed and suggestions for further experimental analyses indicated [7J. 

Up to now all the experimental analyses of B — > decays have been made by heavily 
relying on heavy quark effective theory (HQET) [8] . While HQET provides us with an extremely 
useful tool in understanding and simplifying the non-perturbative dynamics of QCD in the pro- 
cesses involving heavy-light mesons, at the level of precision aimed at the i?-factories, the HQET 
description of the B — > transition matrix element is not as helpful anymore and one should 
attempt doing the QCD analysis instead. 

In the HQET description of this decay, after taking both meson masses to infinity (or m c ^ — > 
oo), the vector form factor F + (q 2 ) -by virtue of the heavy quark flavor symmetry- is related to 
the elastic vector form factor and therefore normalized to one at the zero-recoil. E] To relate that 
normalization to the measured branching fraction one needs to make hazardous computation to 
match HQET with QCD, and include the power corrections that might be uncomfortably large, 
especially those 0(l/m"). Worse even, the symmetry point at which the normalization of the 
form factor is fixed to unity is (7 max = (m B — vtld) 2 [or equivalently w = p B ■ po / l {jnBm B ) = 1], 
where there is no phase space, c l + o^max) = 0; an d therefore the assumptions on the shape of the 
form factor become essential. At first it was believed that the slope of the form factor was enough, 
but later it became clear that an information about its curvature was indispensable [9, 10J. The 
importance of that issue is obvious since the available phase space rapidly grows with w (for lower 
q 2, s). Clearly, a description of this decay that does not rely on HQET is welcome. That statement 
should not be viewed as if HQET is not useful any more. It is still the most valuable framework 
for understanding the dynamics of heavy-light mesons in B — > D** semileptonic and non-leptonic 
decays, and in many other processes, but it is not practical for the exclusive B —> semileptonic 



The term vector (scalar) form factor for F+(q ) (Fo(q )) is related to the fact that in the t-channel it couples 
to the states with quantum numbers J p — 1~ (0 + ). 
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Figure 1: Phase space available from B — > Dtv t is shown in the left plot. Dashed curve corresponds to Co(g 2 ), 
while the solid curve corresponds to c+(q 2 ). The expressions for c5 i +(g 2 ) are given in eq. ([Zj). In the right plot the 
same c+ i0 (q 2 ) are plotted together with c^(q 2 ) (red curve), the phase space available in the case of massless muon 
in the final state. 



modes that are likely to lead to the precision determination of \V c b\, i.e. a determination that 
requires the least number of assumptions about the underlying QCD dynamics. 



2. B — > D£v with minimal theory input: Let us return to eq. ^ and note that the range of 



q s available from this decay is large: 



m 2 < Q 2 < 



[m B - rriD? = 11-63 GeV 2 



(5) 



The form factors FQ + (q 2 ) can be computed on the lattice. The strategy that requires minimum 
assumptions and allows the precision determination exists and it has been implemented in the 
quenched approximation (iVf = 0) |^] The unquenched results with Nf = 2 + 1 dynamical 
light flavors have been recently reported too [13] . An important constraint on the form factors is 
that they are equal at q 2 = 0, 



F+(0) = F (0) 



(6) 



and everywhere else the vector form factor is larger than the scalar one, |i ? +(g ,2 )| > |-Fo(9 2 )|- 
This constraint is useful because the scalar form factor is enhanced by cj(g 2 ), c.f. fig. [TJ and the 
contributions to the decay rate coming from the vector and scalar form factors are competitive 
in size. After taking m 2 = 0, we can write 



mz 



cW) = cl(q 2 ) + Ac + (q 2 ) = c^q 2 ) - A 3 / 2 (g 2 , m%, m 2 D ) ^ 



in: 



(7) 



4 Another strategy has been proposed and implemented recently in the computation of the B-meson decay 
constants /b but not in the computation of B — > D transition form factors [12) - 
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where g^ax = ( m -B ~~ nifl) 2 . On the other hand 
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so that the ratio from eq. ([!]) can be written as 
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Most of these integrals can be evaluated by using the experimentally determined form factor 
|Vd>|Cr(u;), extracted from the differential branching fraction in ref. [14|, and related to \ V c h\Fj r (q 2 ) 
via, 



2 2 . 2 o rp , 2\ ™> B + m D . 

q = m B + m D - 2m B m D w , F + (q ) = G(w) 



w(q 2 ) 



(12) 



For the numerator in the lowest three q 2 bins, each containing a large fraction of events, lead 
to an accurate result, 
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(13) 



For the denominators in i?^ )T we need the form factor F + (q 2 ) for q 2 G [m 2 , which is difficult 

to extract from experiment alone because of the smallness of phase space at larger q 2, s (smaller 
it?), as we can see from fig. [TJ One can instead combine the experimental results obtained within 
q 2 G [m 2 , 8 GeV 2 ], with the lattice QCD results obtained for q 2 G [8 GeV 2 ,^^] [TT1IT3]. 
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which, after using \V c b\ = 0.0411(16) [18] to multiply the lattice results, leads to 

/**?max 

\V cb \ 2 / c^ + {q 2 )\F + {q 2 )\ 2 dq 2 = (26.3 ± 1.0) + (4.9 ± 0.4) = 31.1 ± 1.1 . (15) 

J m 2 

Notice that due to the phase space suppression, the last range of q 2 's in which we used the lattice 
data, makes less than 10% (20%) with respect to the full range of q 2 accessible from this decay 
with muon (r) in the final state. For the numerator in R T the ratio between the scalar and vector 
form factors is needed too. Lattice QCD results of refs. [T3] can be converted to the form 
factors employed here by 



„ f 2n m B +m D m B -m D 
F+{q ) = . h + {w) - h-(w) 



w(q 2 ) 



2\ r I W + 1 , , N W — 1 



Fo(q ) = ^m B m D ■ h+(w) h-(w) ) , (16) 



tub + run niB — tud 



w(q 2 ) 



and after combining them, we see that the ratio Fa(q )/i^(g 2 ) exhibits a linear q 2 behavior, 
which with the intercept fixed by Fq(0) / F+(0) = 1 allows for an accurate determination of the 
slope a, 

F^) =1 - aq > (17) 

as shown in fig. [2] r] We obtain a = 0.020(1) GeV -2 from the values reported in [TT], and a = 
0.022(1) GeV -2 from the results of ref. |13j . i.e. slightly smaller than the naive pole dominance 
model would suggest a = l/m 2 B , = 0.025 GeV -2 , and very close to the result of the model of 
ref. [IS], a = 0.022 GeV" 2 , or a = 0.021(2) GeV" 2 , as obtained in the recent QCD sum rule 



analyses [TBI HZ] • By using a = 0.021(1), and with the help of eq. (17), for the numerator of R T 
we obtain — 12.9 ± 0.7, which finally gives 

R(D) = 0.31 ±0.02, (18) 

which is more than 1- but less than 2-cr below the BaBar result ([!]). 



3. New Physics: As we saw above the value we obtain is consistent with experiment to within 
2a. Since R(D) requires a minimal non-perturbative QCD theory input, it is tempting to check 
the constraints on New Physics that one can infer from the comparison between theory and the 
BaBar result ([!]). In a generic New Physics scenario without right-handed neutrino that preserves 
the lepton flavor universality (LFU), a coupling from this decay to the scalar, vector and tensor 
operators can be described via|^] 

U cS = -V2G F V cb [(c^b)(£ L ^u L ) + g v (ciMhl^ L ) 

+gs(v)(cb)(i R v L ) + g T {n)(ca^b)(i R a» u v L )) + h.c. , (19) 

5 Please note that the plotted values corresponding to the lattice results with N{ = 2 + 1 are obtained from 
h±(w) that we read off from figs. 6 and 7 of ref. [13J . 

6 A coupling to the pseudoscalar and axial operators cannot be studied in this decay as the corresponding matrix 
elements vanish due to parity. 
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Figure 2: In the left panel we show the ratio of B — ¥ D form factors obtained in the lattice QCD simulations 
in quenched approximation (empty symbols) and those in which N{ = 2 + 1 dynamical flavors are included 
(filled symbols) |13] . In the right plot we show the three regions used in eq. ( |11[ ). The same lattice data are used 
in the large q 2 region, while the | V c b\F+(q 2 ) at q 2 < 8 GeV 2 is extracted from the measured differential branching 
fraction for B — >• D\w [14]. 

where the dimensionless couplings gv,s,T ot m %/ w ith mNP being the New Physics scale. 

The differential decay rate from eq. Q now becomes, 
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(D(p')\ca^b\B(p)} = -i {p^ v - v'^Vv) 



2F T (q 2 ,f,) 
mB + mo 



(22) 



The above formulas agree with those reported in ref. [20j. A possibility to discern a small gy 
from this experiment seems very unlikely. It was recently searched while checking the unitarity 
of the first raw of the CKM matrix |19j and was found to be consistent with zero. On the other 
hand gs(n) 7^ is plausible but its value is expected to be very small as the left-right operator 
lifts the helicity suppression and affects both B — > Dtv and B — >■ DfiP decays. Its non-zero 
value could also be a source of difficulties for the D° — D° mixing amplitude, as the left-right 
operators are enhanced by the factor m 2 D /m 2 with respect to the left-left (SM) contribution. A 
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Figure 3: Regions of allowed values for gs(mi,) and grimb), compatible with experimentally measured R(D). 
The first (statistical) error in is treated as Gaussian, while the second (systematic) as uniform. The small region 
within the solid, dashed and dot-dashed white curves correspond to the respective 1-, 2- and 3-<r compatibility with 
both 23 (cxp) (B -> Drv T ) and 23 (cxp) (B -> D\u>y). If the LFU is broken the small region disappears and the above 
contour plot describes gs,T( m b)- The thick dot represents the Standard Model value, namely gs,T(mt) = 0. 



sizable effect could also be seen in D — > Vj decays, as those too are governed by the loops, with 
the down- type quarks propagating in the loop and therefore sensitive to gs(fJ>) 7^ [21] . From 
R(D) alone we get a very loose constraint on gs(fnb) (cf. the contour plot in fig. [3]). Requiring 
the compatibility of the theoretical expression for B(B — > Dfiv^) obtained by using eq. (20) with 
gs( m b) 7^ and gy = gTijnb) = 0, and the measured value [H], restricts the allowed gs(mb) to 
a small region also indicated in fig. [3] For example, if gs(mb) is real then the la compatibility 
with experiment would allow —0.37 < gs( m b) < —0.05, while the 3a compatibility would amount 
to —0.53 < gs(mb) < +0.20. Note that gs(mb), extracted from (tree level) semileptonic process 
should be run to \i = tonp before using it in the loop induced processes. With the help of the 
MS mass anomalous dimension [22], we find gsi^NP = 1 TeV) = 0.58 gsip^b)- In the above 
discussion we assumed a common practice of using the standard quark masses renormalized in 
the MS scheme at \i = nib, the values of which can be found in ref. [23] . 

If in eq. (20) we set gy = gsijnb) = and allow for flr(wib) 7^ 0, then the possible values for the 
real and imaginary parts that are compatible with R{D) [1] are those in the contour plot shown 
in fig. [3} To do that we obviously needed the tensor form factor which has not been computed 
on the lattice nor in the QCD sum rules. To our knowledge it was only computed in the model 
of ref. [15| from which we learn that Fxiq 2 )/ F + (q 2 ) = 1.03(1) is a constant, in agreement with 
naive expectations based on the pole dominance. Again, R(D) alone is not constraining much 
the possible values of grimb), and the compatibility with the measured B{B — > D^v^ [H] helps 
selecting a smaller region, also shown in fig. [3j If Im gTi^b) = 0, we obtain 0.3 < gT^b) ^ 1-5 
and —0.6 < gr^b) < 2.1, from the requirement of respective 1- and 3a compatibility with both 
experimental R(D) and B{B — > Dfiu^). Note again that gxim^p = 1 TeV) = 0.82 grimb), where 
we used the QCD anomalous dimension of the tensor density operator [24j . 
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As we see from the above discussion, one could easily bridge the gap between experimental 
result and the SM prediction for B{B — > Dtu t ) by invoking the New Physics effects that preserve 
LFU. If we give up LFU, the couplings gv,s,T(^) in the Lagrangian |l9| become gysri^)' depen- 
dent on the lepton species too. Many concrete New Physics models, such as various types of Two 
Higgs Doublet Model (2HDM), break LFU and, for example, the scalar coupling g s (/x) becomes 
proportional to the mass of the charged lepton [3]. If so, the constraint from B(B — > D\iv^) is 
a factor of m^/m T less sensitive to New Physics than in the case of LFU illustrated above. The 
net effect would be that the small regions within the white curves depicted in fig. [3] would simply 
disappear, and the preferred regions of the scalar coupling g T s {mi,) would be those represented by 
the blue contours in the left plot of fig. [3j 



4. Concluding comments: In closing this letter we would like to make the following comments: 

• Assuming the lepton flavor universality, which has been experimentally verified to a very 
good accuracy [23] , the result ([I]) could be an indication of New Physics, if incompatibility 
with the Standard Model is indeed shown to be significant (more than 3cr). That test 
[of compatibility with the Standard Model] can be made experimentally, with a minimal 
theory input, as discussed in this letter. Here we used the lattice QCD results for F + {q 2 ) at 



larger q 2, s because the full branching fractions were reported in ref. [JJ. Our value (18) is 
compatible with experiment within less than 2a. If, instead of measuring the full branching 
fractions for both decay modes, the experimenters made a cut at about q 2 « 8 GeV 2 , then 
the full shape of the needed vector form factor could be reconstructed from the differential 
branching fraction of B — > D\xv j3j H] . The only information needed from theory is then the 



ratio of the scalar and vector form factors (17), which is quite accurately known from lattice 



QCD, with the values that agree with quark models, and with recent QCD sum rule studies. 
We hope such an analysis will be made by both BaBar and Belle. By using the vector form 
factor data from ref. [14] only, and by integrating the decay rates up to q 2 ut = 8 GeV 2 , we 
obtain 



B(B -> Dtv t ) 



B(B -> D^) 



= 0.20 ± 0.02 . (23) 

g 2 <8 GeV 2 



• B — > D^jbU decay can be the mode allowing the most reliable extraction of \ V c b\ as it requires 
the least number of assumptions. A discussion made in ref. |14j showed that one can find 
a range of q 2, s in which both the experimental and lattice QCD errors can be kept small 
and therefore allow for a very clean extraction of |V^,|. As for the large measured value 
for B(B — > Dtu t ), an independent measurement by Belle is necessary. A measure of its 
partial branching fraction would also be helpful to permit a direct comparison with the 
lattice QCD results for the form factors. In both CclSGS, Si special care should be devoted to 
the systematics related to the presence of B — > D£i"y events in the selected sample, with 
photon being soft. Such events affect the neutral and the charged i3-meson semileptonic 
decays differently [2]. 

• As shown above, the measured R{D) and B(B — > Dfiu^) provide quite good constraints on 
the new physics couplings gs,T(mb), that then can be used at the loop induced processes, 
after running to 55,t(^np)- 
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• An experimental study of B s — > D s iv decay rates would be even more advantageous. It 
would eliminate the chiral extrapolations of the lattice results in the valence light quark 
mass. A fully unquenched lattice QCD study along the lines presented in ref. [H] would be 
very welcome too. It would be an important independent check of the results of ref. |13| . 
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